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Y^' Abstract 

D , So far no theoretical tool for the comprehensive description of exclusive electropro- 

i-^ ' duction of vector mesons off nuclei at medium energies has been developed. We suggest a 

light-cone QCD formalism which is valid at any energy and incorporates formation effects 
K> , (color transparency), the coherence length and the gluon shadowing. At medium energies 

5-H I color transparency (CT) and the onset of coherence length (CL) effects are not easily 

separated. Indeed, although nuclear transparency measured by the HERMES experiment 
rises with Q^, it agrees with predictions of the vector dominance model (VDM) without 
any CT effects. Our new results and observations are: (i) the good agreement with the 
VDM found earlier is accidental and related to the specific correlation between Q^ and 
CL for HERMES kinematics; (ii) CT effects are much larger than have been estimated 
earlier within the two channel approximation. They are even stronger at low than at high 
energies and can be easily identified by HERMES or at JLab; (iii) gluon shadowing which 
is important at high energies is calculated and included; (iv) our parameter-free calcula- 
tions explain well available data for variation of nuclear transparency with virtuality and 
energy of the photon; (v) predictions for electroproduction of p and (p are provided for 
future measurements at HERMES and JLab. 



1 Introduction: interplay of absorption and shadowing 

1.1 Color transparency 

The nuclear medium is more transparent for colorless hadronic wave packets than predicted 
by the Glauber model. One can treat this phenomenon either in the hadronic basis as a 
results of Gribov's inelastic corrections ||l|, or in QCD as a result of color screening P, Q, an 
effect called color transparency (CT) (see also the review 0]). Although the two approaches 
are complementary, the latter interpretation is more intuitive and straightforward. Indeed, a 
point-like colorless object cannot interact with external color fields, therefore its cross section 
vanishes a{r) oc r^ at r — i> [§]. This fact naturally explains the correlation between the cross 
sections of hadrons and their sizes P, |^, 0| . When a colorless wave packet propagates through 
a nucleus, the fluctuations with small size have an enhanced survival probability which leads 
to a non-exponential attenuation oc 1/L [0], where L is the path length in nuclear matter. 

Diffractive electroproduction of vector mesons off nuclei is affected by shadowing and ab- 
sorption which are different phenomena. Final state absorption of the produced meson exists 
even in the classical probabilistic approach which relates nuclear suppression to the survival 
probability W{z, h) of the vector meson produced at the point with longitudinal coordinate z 
and impact parameter 6, 



W{z, h) = exp 
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where pA^b, z) is the nuclear density and aJJ^ is the inelastic VN cross section. 

Going beyond the VDM one realizes that the diffractive process initiating the production 

of the vector meson on a bound nucleon is 7* A^ -^ qq N (with possible glue). A photon of 

high virtuality Q^ is expected to produce a pair with a small ~ l/Q"^ transverse separation^. 

The basic idea of CT is that such small size should lead to a vanishing absorption when the 

colorless qq wave packet propagates through the nucleus. However, the pair may evolve in size 
^In fact, the situation is somewhat more com.phcated. For very asymmetric pairs when the q or q carry 
almost the whole photon momentmn, the pair can have a large separation, see Sect. |2.2| 



during the time of propagation due to transverse motion of the quarks. Besides, the medium 
filters out large size configurations which have larger absorption cross section|, an effect known 
as color filtering. Eventually, the resulting distribution amplitude of the qq wave packet must 
be projected onto the wave function of the V meson. 

The time scale characterizing the evolution of the qq wave packet can be estimated based 
on the uncertainty principle. One cannot decide whether the ground state V is produced or 
the next excited state V, unless the process lasts longer than the inverse mass difference. In 
the rest frame of the nucleus this formation time is Lorentz dilated, 

tf = -r-2 2 ' (2) 

my — rriy 
where z/ is the photon energy. 

A rigorous quantum-mechanical description of the pair evolution was suggested in |^ and is 
based on the light-cone Green function technique. This approach is presented below in Sect. Q. 

A complementary description of the same process in the hadronic basis looks quite different 
ip. The incident photon may produce different states on a bound nucleon, the V meson ground 
state or an excited state. Those states propagate through the nucleus experiencing multiple 
diagonal and off-diagonal diffractive interactions, and eventually the ground state is detected. 
According to quark-hadron duality we expect these two descriptions to be equivalent. In prac- 
tice, however, neither of them can be calculated exactly, and therefore each has advantages and 
shortcomings. For example, electroproduction of light vector mesons on a nucleon cannot be 
calculated perturbatively without reservations, while in the hadronic basis one can make use 
of experimental data which include all nonperturbative effects. On the other hand, for excited 
meson states no data are available for the diagonal and off-diagonal diffractive amplitudes which 
one can estimate in the quark representation. The two approaches are complementary, they 
rely on different approximations and their comparison may provide a scale for the theoretical 

uncertainty involved. 

^Absorption does not mean disappearance or stopping of the quarks. High energy partons usuaUy lose only 
a very small (energy independent) fraction of their energy, primarily via soft QCD processes. Absorption means 
color-exchange interaction which switches from the exclusive channel to an inclusive process. 



1.2 Effects of coherence: shadowing of quarks and gluons 

Another phenomenon, shadowing, is also known to cause nuclear suppression. In contrast to 
final state absorption, it is a pure quantum-mechanical effect which results from destructive 
interference of the amplitudes for which the interaction takes place on different bound nucleons. 
It can be interpreted as a competition between the different nucleons participating in the 
reaction: since the total probability cannot exceed one, each participating nucleon diminishes 
the chances of others to contribute to the process. 

The cross section of photoproduction is very small since it includes the fine structure con- 
stant. Applying the Glauber formula one should expect no visible shadowing. However, this 
is true only at low energies. It has been realized back in the 60s (see the review [|l^) that the 
photon interacts via its hadronic fiuctuations. Therefore, if a fiuctuation can propagate over a 
distance comparable or longer than the nuclear radius, it may interact with a large hadronic 
cross section which causes shadowing. The small probability to create such a fiuctuation enters 
only once, otherwise the fiuctuation interacts strongly. Thus, the fiuctuation lifetime provides 
the time scale which controls shadowing. Again, it can be estimated relying on the uncertainty 
principle and Lorentz time dilation as, 

_ 2 z/ 

^^ - g^T< • ^^^ 

It is usually called coherence time, but we also will use the term coherence length (CL), since 
light-cone kinematics is assumed, Ic = tc (similarly, for formation length If = tf). CL is 
related to the longitudinal momentum transfer q^ = l/lc in 'y* N ^ V N, which controls the 
interference of the production amplitudes from different nucleons. 

Initial state shadowing indeed has been observed in many reactions where no final state 
absorption is expected, for example in the total photoabsorption cross section on nuclei (see 
|T0[), the inclusive deep-inelastic cross section ||TT], |12|, the total neutrino- nucleus cross section 



||T3| , the Drell-Yan reaction of dilepton production [|1J, |T5| , etc. In the case of electroproduction 
of vector mesons off nuclei shadowing and absorption happen with the same cross section which 
makes it difficult to disentangle the two sources of nuclear suppression. Nevertheless, it is easy 



to identify the difference in the two hmiting cases which we illustrate for the example of the 
vector dominance model (VDM). The first case is the limit of small Ic, shorter than the mean 
internucleon spacing ~ 2 fm. In this case only final state absorption matters. The ratio of the 
quasielastic (or incoherent) 'y* A —* V X and 'y* N -^ V X cross sections, usually called nuclear 
transparency, reads, 
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In the limit of long l^ it takes a different form. 



y^mc 



Ic'»Ra 



(5) 
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where we assume a^i^ <^ o"]^^ for the sake of simplicity. TA{b) is the nuclear thickness function 

oo 

TA{b)= J dzpA{b,z) . (6) 
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The exact expression beyond VDM which interpolates between the two regimes (^ and (|^) can 



be found in [16 



One can see that the V meson attenuates along the whole nucleus thickness in Eq. (^), 
but only along roughly half of that length in (^. This confirms our conjecture that nuclear 
shadowing also contributes to (^ increasing suppression. This may be also interpreted as an 
analog of the quark nuclear shadowing measured in DIS off nuclei, but the absorption effects 
make this analogy rather shaky. 

Gluon shadowing also suppresses electroproduction of V mesons. Different (but equivalent) 
descriptions of gluon shadowing are known. In the infinite momentum frame of the nucleus 
it looks like fusion of gluons which overlap in longitudinal direction at small x, leading to a 
reduction of gluon density. In the rest frame of the nucleus the same phenomenon looks as a 
specific part of Gribov's inelastic corrections IQ. The lowest order inelastic correction related to 
diffractive dissociation V N ^ X N [ITTtl contains PPR and PPP contributions (in terms of the 



triple- Regge phenomenology, see [1^)- The former is related to the quark shadowing already 
discussed above, while the latter, the triple-Pomeron term, corresponds to gluon shadowing. 
Indeed, only diffractive gluon radiation can provide the Mx dependence dadd/dMx oc 1/M| of 
the diffractive dissociation cross section. 

In terms of the light-cone QCD approach the same process is related to the inclusion of 
higher Fock components, \qqnG), containing gluons [1^. Such fluctuations might be quite 
heavy compared to the simplest \qq) fluctuation, therefore, they have a shorter lifetime [^ and 
need a higher energy to be relevant. 

1.3 Outline of the paper 

In Sect. ^ we present the light-cone (LC) approach to diffractive electroproduction of vector 
mesons in the rest frame of the nucleon target. The central issue of this approach, the universal 
interaction cross section for a colorless quark-antiquark dipole and a nucleon, is presented in 



Sect. p.l| . It cannot be reliably evaluated theoretically and is fitted to the data for the proton 
structure function in a wide range of xbj and Q^. 

The LC wave function for a quark-antiquark fluctuation of the virtual photon is presented 
in Sect. |2.2| for both, free and interacting qq pairs. In the latter case we apply the LC Green 
function approach and introduce into the two-dimensional Schrodinger equation a nonpertur- 
bative real LC potential describing the qq interaction. The model for the LC wave of a vector 



meson is described in Sect. 2.3. 



As a rigorous test of the model we calculate in Sect. |2^ the cross section of elastic electro- 
production of p and mesons off a nucleon target. These parameter- free calculations reproduce 
both energy and Q^ dependence remarkably well, including the absolute normalization. Since 
we use the nonperturbative LC photon wave function it is legitimate to do calculations down 
to Q^ = 0. Agreement with data for real photoproduction of p and (p is also good. 

Sect. ^ is devoted to incoherent production of vector mesons off nuclei. In Sect. p.l| the 
Green function describing propagation of a qq in the nuclear medium is modified to incorporate 



absorption. This is done by introducing an imaginary part of the potential into the two- 
dimensional LC Schrodinger equation for the Green function. Different limiting cases of short 
and long coherence and formation lengths are considered. The central results of the paper 
is Eq. ( p2D for the cross section of incoherent vector meson production in the most general 



case. Numerical calculations and comparison with available data are presented in Sect. ^ 
Nuclear transparency turns out to be a result of a complicated interplay between coherence 
and formation length effects. Although variation of Ic with Q^ can mimic CT at medium and 
low energies, one can map experimental events in Q^ and z/ in such a way as to keep Ic = const. 
Unexpectedly, the exact solution found in the present work is very different from the two- 
coupled-channel approximation of f^ and predicts a much more pronounced effect of CT. This 
makes it feasible to find a clear signal of CT effects in exclusive production of p mesons in the 
current and planned experiments at HERMES and JLab. 

Coherent production of vector mesons off nuclei leaving the nucleus intact is studied in 
Sect. ^. The formalism described in Sect. |4.1| is simpler than in the case of incoherent pro- 



duction. The detailed calculations and the comparison with data are presented in Sect. O 
The effect of CT on the Q^ dependence of nuclear transparency at Ic = const is weaker than 
in the case of incoherent production and is difficult to be detected at low energies since the 
cross section is small. Our results for the differential cross section of coherent production of p 
demonstrate also a weak sensitivity to the CT effects. 

Besides CL, there are other effects considered in Sect. ^ which can mimic the phenomenon of 
CT. First, the standard lowest order inelastic corrections well fixed by available data are known 
to make the nuclear medium more transparent at higher energies. Since i/ is a rising function 
of Q^ at fixed Ic, nuclear transparency increases with Q^. These corrections are estimated in 
Sect. ^TT| and the effect is found to be too weak to mock CT. Another source of rising Q^ 
dependence of the nuclear transparency is the finite p lifetime which might be important at low 



energies. This effect evaluated in Sect. |5.2| is also found to be negligibly small. 

Exclusive production of vector mesons at high energies is controlled by the small-x^j physics, 
and gluon shadowing becomes an important phenomenon. It affects the cross section of inco- 



herent vector meson production in a two-fold way. While the production of V^ on a bound 
nucleon is suppressed, the nuclear medium becomes more transparent enhancing the survival 
probability of qq wave packets traveling through the nucleus. At the same time, the cross sec- 
tion of coherent production can be only diminished. In Sect. ^ gluon shadowing is calculated 
and included in the calculations for nuclear transparency. 

The results of the paper are summarized and discussed in Sect. 0. An optimistic prognosis 
for the CT discovery potential of future experiments at HERMES and JLab is made. 

2 Light-cone dipole phenomenology for elastic photopro- 
duction of vector mesons 7*Ar — > V N 

In the light-cone dipole approach the amplitude of a diffractive process is treated as elastic 
scattering of a qq fluctuation of the incident particle. The elastic amplitude given by convolution 
of the universal flavor independent dipole cross section for the qq interaction with a nucleon, 
(jqq, which is introduced in |^, and the initial and final wave functions P|. Thus, the forward 
production amplitude for the exclusive photo- or electroproduction of vector mesons 7*A^ — »• 
V N can be represented in the form 

1 
M,,m^vn{s,Q'') = {V\ag{p,sM) = jdaj(fr vl/;^(f, a) a^-.lf, s) vi/^*(f, a, Q^) (7) 



with the normalization 

da 
It 



t=0 16 TT 



l-^l'^ (8) 



In order to calculate the photoproduction amplitude one needs to know the following in- 
gredients of Eq. (0): (i) the dipole cross section agg{r,s) which depends on the qq transverse 
separation r and the cm. energy squared s. (ii) The light-cone (LC) wave function of the 
photon \E'^. (r*, a, Q^) which also depends on the photon virtuality Q^ and the relative share a 
of the photon momentum carried by the quark, (iii) The LC wave function \E'y(r, a) of the 
vector meson. They are presented in the following sections. 
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Note that in the LC formahsin the photon and meson wave functions contain also higher 
Fock states |gg), \qqG), \qq2G), etc. Should one add their contribution to Eq. (|])? A word of 
caution is in order. The energy dependence of the total cross section including the dipole one 
as is given in Eq. (J^), originates from inelastic collisions with gluon bremsstrahlung, a process 
related to the forward elastic amplitude via unitarity. Those inelastic collisions also can be 
described in terms of the Fock components containing gluons. Thus, one would double count 
if both the energy dependent dipole cross section and the higher Fock states were included. 
Either one should rely upon the Fock state decomposition treating interaction of each of them 
as energy independent, or one should restrict ones consideration to the lowest \qq) component, 
but implicitly incorporate the effects of higher Fock states into the energy dependence of the 
dipole cross section aqq. We stand with the latter approach in the present paper. 

However, as for nuclear targets, one must explicitly include into ones consideration the 
higher Fock states because their eikonalization leads to gluon shadowing. We come back to this 
problem in Sect. ^ 

2.1 Phenomenological dipole cross section 

The cross section o"gg(r*, s) for the interaction of a qq dipole of transverse separation r with 
a nucleon, first introduced in [0 , is a flavor independent universal function of r and energy. 
It allows to describe in a uniform way various high energy processes. This cross section still 
cannot be predicted reliably because of poorly known higher order pQCD corrections and 
nonperturbative effects. However, it is known to vanish quadratically o"gg(r, s) oc r^ as r — > 
due to color screening, a property usually called color transparency. On the other hand, one 
may expect the dipole cross section to level off at large separations. This may happen if the 
quark density in the proton already saturates in the xbj range of HERA pT| , p2| . Whether 
this was already observed at HERA in the small xbj domain is difficult to say. One can fit 



the data perfectly either assuming saturation |22|, or with the pure DGLAP evolution. One 



can interpret the leveling off of the dipole cross section not only in terms of saturated parton 



density (these two are not identical). Another scenario relates the flat behavior of o"qg(r, s) at 



large r to the averaged gluon propagation length tq. For r^ ^ Vq one arrives in this case at the 
additive quark model: the dipole cross section is a sum of quark-quark cross sections, i.e. the 
<Jqq{r) levels off at large separations. 

At small separations the dipole cross section should be a function of r and xbj ~ l/(r^s) 
to reproduce Bjorken scaling. A corresponding simple and popular parameterization has been 
suggested in [^]. It well describes data for DIS at small x and medium and high Q^. However, 
at small Q^ it cannot be correct since it predicts energy independent hadronic cross sections. 
Besides, xbj is not any more a proper variable at small Q^ and should be replaced by energy. 
Since we want our approach to be valid down to the limit of real photoproduction, we choose 
the parametrization suggested in |23| which is similar to one in ||22|, but contains an explicit 
dependence on energy. 



(^gqir,s) 



(^o[s) 



I _ g-'-V^oCs) 



(9) 



It correctly reproduces the hadronic cross sections for the choice 



o-o(s) = aZ{s) 



1 + 



-0.14 



mb 



(10) 



ro(s) =0.88(— j f m . (11) 

Here {rlf^) = 0.44 fm^ is the mean pion charge radius squared; Sq = 1000 GeV^. The cross 
section (Ttoti^) was fitted to data in 123, |2H|, 

/ „ \0.079 / „ \ -0.45 

ar.^(.) = 23.6 - +0.032 - mb . (12) 

VSq/ \So/ 

It represents the Pomeron and Reggeon parts corresponding to exchange in of gluons and qq, 
respectively. Only the former has been used in the dipole cross section Eq. (|^) to fit the data 
for the proton structure function at small xbj- Unfortunately, the Reggeon part of the dipole 
cross section is poorly known. To the best of our knowledge no phenomenology of it has been 
developed so far. The energy dependence of the Reggeon dipole cross section at small Q^ (or 
x-dependence at high Q"^) dictated by Regge phenomenology is approximately oc 1/^/s (oc y/x). 
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Thus, we can expect Q^ independence of the the exponent in the second term on the r.h.s of 
Eq. (0). The r dependence of this term in the dipole cross section is less known. In order 
to reproduce Bjorken scahng we must assume that the Regge term vanishes at r ^ in the 
same way as the Pomeron part of the cross section, aR{r, s) oc r^. For the sake of simphcity we 
therefore assume the same parameterization as is used for the Pomeron part, Eq. (|^). Then, 
one can just incorporate the Reggeon term into cro(s) as is done in Eq. (|T^). 

Incorporating Reggeons into the LC dipole formalism for nuclear shadowing one should be 



careful with the treatment of multiple interaction terms which have a nonplanar nature \2t\. The 
Reggeon exchange described by planar graphs should not participate in the multiple-scattering 
expansion of the eikonal exponential. 

The dipole cross section Eqs.(||) - ( |T2|) provides the imaginary part of the elastic amplitude. 
It is known, however, that the energy dependence of the total cross section generates also a 



real part p? 



a5,(r, s)^U-i^ ah^ ) ''^''^''' ^^ ^^^^ 

The energy dependence of the dipole cross section Eq. @ is rather steep at small r leading to a 
large real part which should not be neglected. For instance, the photoproduction amplitude of 
7A^ -^ J/'^N rises oc s*^'^ and the real-to-imaginary part ratio is over 30%. At medium energies 
also the Reggeon contribution to electroproduction of light mesons contributes to the real part 
of the elastic amplitude. The replacement Eq. (p!3| ) takes care of it as well, and we use this 
form everywhere in what follows, unless specified otherwise. 

Note that the improvement compared to p2| at large separations leads to a worse description 
of the short-distance part of the dipole cross section which is responsible for the behavior of 
the proton structure function at large Q^. To satisfy Bjorken scaling the dipole cross section 
at small r must be a function of the product sr which is not the case for the parametrization 
in Eq. @. Indeed, the form of Eq. @ successfully describes data for DIS at small x only up 
to Q^ ~ 10 GeV^, and does a poor job at larger values of Q^. Nevertheless, this interval of Q^ 
is sufficient for the purpose of the present paper which is focused on production of light vector 
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mesons at small and moderate Q"^ < 10 GeV^. 

2.2 The qq wave function of the photon 

The perturbative distribution amplitude ( "wave function" ) of the qq Fock component of the 



photon is well known [^, ^, ^], and for transversely (T) and longitudinally (L) polarized 
photons it has the form, 



*^^"(r» = ^^^iF^ Z,xd^^^xK,{er) (14) 



where x ci-nd x are the spinors of the quark and antiquark, respectively; Zq is the quark charge; 
Nc = 3 is the number of colors. KQ{e.r) is a modified Bessel function with 

e^ = a{l-a)Q'^ + ml, (15) 

where mq is the quark mass, and a is the fraction of the LC momentum of the photon carried 
by the quark. The operators O^'^ read, 

O^ = mg (J ■ e + z (1 - 2a) {a ■ n) (e ■ V,-) + (a x e) • V,- , (16) 

O^ = 2Qa{l-a){a-n) . (17) 

Here V,. acts on transverse coordinate r; e is the polarization vector of the photon and n is a 
unit vector parallel to the photon momentum. 

The transverse qq separation is controlled by the distribution amplitude Eq. ([T^ with the 

mean value, 

^^ _ = 18 

e ^Q2^(l-a) + m2 

In pQCD the quarks are treated as free, and one may wonder why they do not fly apart but 
form a wave packet of finite size. It is interference of qq waves produced at different points 
which keeps the transverse separation finite. To reach a large separation the qq pair must 
be produced sufficiently long in advance, longer than the coherence time Eq. (^, such that 
fluctuations loose coherence. Treating the coherence time as lifetime of the fluctuation, one can 
also say that the fluctuation does not have enough time to fly apart. 
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For very asymmetric qq pairs with a or (1 — a) < m^/Q'^ the mean transverse separation 
(r) ~ l/rriq becomes huge since one must use current quark masses within pQCD. A pop- 
ular recipe to fix this problem is to introduce an effective quark mass rrieff ~ ^qcd which 
should represent the nonperturbative interaction effects between q and q. It is more consis- 
tent, however, and straightforward to introduce this interaction explicitly. The corresponding 
phenomenology based on the light-cone Green function approach has been developed in p3| . 

The Green function Gqq{zi,ri; Z2,r2) describes the propagation of an interacting qq pair 
between points with longitudinal coordinates Zi and Z2 and with initial and final separations 
r*! and r2. This Green function satisfies the two-dimensional Schrodinger equation. 



i-j—Ggq{zi,ri]Z2,r2) 

dZ2 



+ Vqq{z2,r2,a) 



Gqq{zi,ri;Z2,r2) 



(19) 



2 z/a (1 — a) 

Here u is the photon energy. The Laplacian A^ acts on the coordinate r. 

The imaginary part of the LC potential Vqq{z2,r2,C() in (0) is responsible for attenuation 
of the qq in the medium, while the real part represents the interaction between the q and q. 
This potential is supposed to provide the correct LC wave functions of vector mesons. For the 
sake of simplicity we use the oscillator form of the potential, 

a^(a)f2^ 



ReVqq(z2,r2,a) 



(20) 



2 i^a(l — a) 

which leads to a Gaussian r-dependence of the LC wave function of the meson ground state. 
The shape of the function a{a) will be discussed in the next section. 

In this case equation (|19D has an analytical solution, the harmonic oscillator Green function 



0, 



Gqq(zi,fi; Z2,r2 



a^(a) 



2 Tx i sinfu; Az 



■exp 



ia?{a.) 
smiio Az) 



{r\ + rl) cos(c<j A2;) — 2 n • r2 
ie 



xexp 



'2A2 



2z/a (1 — a) 



where Az = Z2 — zi and 



uo 



a?{a) 
V a{\ — a) 



(21) 



(22) 



13 



The boundary condition is Ggq{zi,ri; Z2,r2)\z2=zi = ^^(^i — "^2). 

The probabihty amphtude to find the qq fiuctuation of a photon at the point Z2 with 
separation f is given by an integral over the point zi where the qq is created by the photon 
with initial separation zero, 



*«'(*'• "' = i.'E^'a) J "'' (^- '^"■'*) '''-''^" *''• ''-' '> 



ri=0 



(23) 



The operators O^'^ are defined in Eqs. (|TB|) and (0). Here they act on the coordinate fi. 
If we write the transverse part as 



xO^X = A + B-\/r, 



(24) 



then the distribution functions read, 



m(f,a) = Zg,/d^ A$o(e,r,A) + S$i(e,r,A) 



*^ (f, a) = 2 Zg^/a^Q a{l - a) x ^ ■ n x %{^, r, A) 



where 



A 



2aHa) 



The functions $0,1 in Eqs. ( p5|) and (|26|) are defined as 

AeV2 



*°^^'^'^) = i/^';hM"^ 



e^r 



<l>i(e,r,A) = --/rft 
on J 




A 



sh(At) 



exp 



4 
AeV2 



cth(At) - 1 
cth(At) - 1 



(25) 
(26) 

(27) 

(28) 
(29) 



Note that the q — q interaction enters Eqs. ( 125|) and (^) via the parameter A defined in 
|27|) . In the limit of vanishing interaction A ^ (i.e. Q^ -^ 00, a is fixed, a 7^ or 1) Eqs. (p5|) 



produce the perturbative expressions of Eq. (p^). 
With the choice a^(a) oc a{l ~ a) the end-point behavior of the mean square interquark 
separation (r^) oc l/a{l — a) contradicts the idea of confinement. Following [^ we fix this 
problem via a simple modification of the LC potential. 



a^{a) = Oq + 4a^ a(l — a) 



(30) 
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The parameters oq and Oi were adjusted in [23| to data on total photoabsorption cross section 
32| , |33[| , diffractive photon dissociation and shadowing in nuclear photoabsorption reaction. 



The results of our calculations vary within only 1% when oq and ai satisfy the relation, 



,1.15 



(0.112)2 Qgy2 



'1-v 



1.15 



(0.165)2 GeV^ 



(31) 



where v takes any value < i; < 1. In view of this insensitivity of the observables we fix the 
parameters at f = 1/2. We checked that this choice does not affect our results beyond a few 
percent uncertainty. 



2.3 The meson wave function 

To describe electroproduction reactions it is natural to work in the infinite momentum frame 
of the virtual photon and use the LC variables for the qq pair, the transverse separation f and 
the fraction a = pt /py of the total LC momentum carried by the quark. The wave functions 
of light vector mesons are poorly known both in the rest and infinite momentum frames. A 



popular prescription [Q is to apply the Lorentz boost to the rest frame wave function assumed 
to be Gaussian which leads to radial parts of transversely and longitudinally polarized mesons 
in the form, 

^ ,. r n/(l — n/\ r^ 

(32) 



$^'^(f,a) = C^'^a(l- 


a)/(a)exp 


— 


2i?2 


with a normalization defined below, and 


/(a) = exp 


r 2 n2 1 




2a{l-c 


M 





(33) 

This procedure is ill motivated since the qq are not classical particles. As a result of the boost to 
the infinite momentum frame many new Fock components are created. Nevertheless, a detailed 
analysis of this problem |3^ leads to the same form (|32|) which we use in what follows with the 
parameters from |36|, R = 0.59 fm and mg = 0.15 GeV. 
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We assume that the distribution amphtude of qq fluctuations for the vector meson and for 
the photon have a similar structure ||36|. Then in analogy to Eqs. ( ^Sf ) - (|26|) , 



^^(r,a) = (A + 5-V)$^(r,a); 



(34) 
(35) 



Correspondingly, the normalization conditions for the transverse and longitudinal vector 
meson wave functions read, 



Nc I d r I daim. 



$^(f,a) + a^ + (l-a) 



dr^Ur,a] 



ANc d^r daa^{l-aymy $y(f,a 



1. 



(36) 
(37) 



2.4 Cross section on a nucleon, comparison with data 

Now we are in the position to calculate the forward production amplitude 7* A^ -^ V N for 
transverse and longitudinal photons and vector mesons using the nonperturbative photon wave 
functions Eqs. (^), (^) and for the vector meson Eqs. (|3^), (|35|) . We verify the LC approach 
by comparing with data for nucleon target. This is a rigorous test since we have no free 
parameters. 

The forward scattering amplitude reads. 



M 



■yN^V N 



:s,Q' 



t=0 



Nc Zg y/a^J d^r aqq{r, s) J da\ml <|)o(e, f, A)^y(f, a) 


$i(e,f,A)-V.*^(f,a)}; (3^ 



+ W + il-af 



•M-^N^V N^ 



s,Q') 



t=0 



X 



4 Nc Zq ^/a^my Q / d'^r aqqif, s) 
1 
/ da c? (1 - of $o(e, f, A)^y(f, a) . 



(39) 



IGvr (icr^' jdt . We include the real part of 



These amplitudes are normalized as |A^"^'^p — t«--'^-- ' 

the amplitude according to the prescription described in Sect. ^]l|. In what follows we calculate 

the cross sections o = o^ ^ ta^ assuming that the photon polarization is e = 1. 
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Now we can check the absolute value of the predicted cross section by comparing with data 
for elastic electroproduction 'j* p ^ V p for p and (p mesons. Unfortunately, data are available 
only for the cross section integrated over t, 



a'^'^{-f*N ^VN) 



(40) 



where the t-slope of the differential cross section which cannot be properly predicted by the 
approach under consideration. Our strategy is to predict the numerator in P0|), and compare 
with data for the cross section and the slope. 

Our predictions are plotted in Fig. [^ together with the data on the Q^ dependence of the 



cross section from NMC, HI and ZEUS P, p8|, p9|, m). 
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Figure 1: Q"^- dependence of the cross section for the reactions j* p ^ pp 
(left) and 'j* p ^ (f)p (right). The dashed and sohd curves are compared 
with data at W = 15 GeV ^ and at 75 GeV (^, |3|/ for p and g^ for 
(j)), respectively. 



We use the Q^ dependent slope of the differential cross section da{Y^ ~^ VN)/dt oc 



exp 



B^j^iQ"^) t parametrized as [^ 



B,.n{s,Q')=P^{s) + 



Q^ + ml 



2 



' Q'^ + ni 



m 



(41) 



V 
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A fit to the data [37, 42 1 from fixed target experiments for the Q^ dependent slope in p pro- 
duction give the parameters for PF ^ 10 - 15 GeV, (3^ = (6.2 ± 0.2) GeV'^, (3{ = 1.5 ± 0.2. 

I, E6|, E^ for p production at W = 75 GeV we get, 



Using data from HERA ^ |3|, ||, ||, 
f3^ = (7.1 ± 0.1) GeV"^ /?f = 2.0 ± 0.1. 

Repeating the same analysis for production we get at VT ~ 10 — 15 GeV from data [p^, HS 



/3o* = (5.9±0.1) GeV-^ /3f = 0.5±0.1. Data from HERA |9|, |g give, (3^ = (6.7±0.2) GeV'^ 
Pf = 1.0 ± 0.1. For calculations shown in Fig. |^ we use the central values of these parameters. 
Our approach which includes the effects of the nonperturbative interaction between the q 
and q in the photon fluctuation is designed to describe the low Q'^ region as well. To test 
it we compare with data [^0|, ^, Q ^ ^ for the energy dependence of the cross section 
of real p photoproduction in Fig. ^. We use the energy dependent slope parameter, B^j^ = 
B^ + 2a'ln(s/so) with a' = 0.25 GeV"^ and B^ = 7.6 GeV'^ Sq = 20 GeV^ fitted to data 
53| , pO| , ^ , ^ ^. The Pomeron part of the dipole cross section depicted by the dashed curve 



in Fig. 1^ cannot explain the data at low energies, W^ < 15 GeV, while the addition of the Regge 
term (solid curve) leads to a good agreement for all energies. We also found a good agreement 
with data for real photoproduction of 0, but skip the comparison since there are very few data 
points. 

The normalization of the cross section and its energy and Q^ dependence are remarkably well 
reproduced in Figs. [^ - 0. This is an important achievement since the absolute normalization 
is usually much more difficult to reproduce than nuclear effects. For instance, the similar, but 
simplified calculations in [^ underestimate the J/^ photoproduction cross section on protons 
by an order of magnitude. 

As a cross-check for the choice of the p° wave function in Eqs. (0) and (^0]) we also 
calculated the total p°-nucleon cross section, which is usually expected to be roughly similar to 
the pion-nucleon cross section a^^ ~ 25 mb. The p-nucleon total cross section has the form. 



a 



pN 
tot 



N, 



c 



cPr / da { mt 



$^(f,a) 



+ \a' + (1 



a] 



dr^{r{r,a) 



(^qq\r,S} 



(42) 



We calculated alot with the p meson wave function in the form Eq. (|32|) with the parameters 
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Figure 2: Energy dependence of the real photoproduction cross section on 
a nucleon, 'j p -^ p^ p. Our results (solid curve) are compared with data 

and ZEUS MB, 



from the fixed target J5C| , |5J|/ , and collider HERA HI 

P^/ experiments. The dashed curve contains only the gluonic exchange 

in the t-channel. 



described in the Sect. 2.3. For the dipole cross section we adopt the KST parameterization 



which is designed to describe low-Q^ data. Then, at z/ = 100 GeV we obtain a 
which is quite a reasonable number. 



pN 
tot 



27 mb 



3 Incoherent production of vector mesons off nuclei 

In diffractive incoherent (quasielastic) production of vector mesons off nuclei, 7* A —* V X, 
one sums over all final states of the target nucleus except those which contain particle (pion) 
creation. The observable usually studied experimentally is nuclear transparency defined as 



Tr'X" 



a. 



-fA^VX 



Aa^* 



'y*N^VN 



(43) 
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The t-slope of the differential quasielastic cross section is the same as on a nucleon target. 
Therefore, instead of integrated cross sections one can also use the forward differential cross 
sections Eq. (||) to write, 

(44) 



y^mc 



1 
A 



j'-A^VX 



:s,Q' 



M^* 



■y*N^VN 



:s,Q') 



3.1 The LC Green function approach 



One should decompose the physical photon I7*) into different Fock states, namely, the bare 
photon |7*)o, Iqi), IqqG), etc. The higher states containing gluons are vital to describe the 
energy dependence of the photoproduction reaction on a nucleon. As far as nuclear effects are 
concerned, those Fock components also lead to gluon shadowing. However, as we mentioned 
above, these fluctuations are heavier and have a shorter coherence time (lifetime) than the 
lowest \qq) state. Therefore, at medium energies only \qq) fluctuations of the photon matter. 
Gluon shadowing related to the higher Fock states will be considered later. 

Propagation of an interacting qq pair in a nuclear medium is described by the Green function 
satisfying the evolution Eq. (|1^). However, the potential in this case acquires an imaginary 
part which represents absorption in the medium (see ([l|) for notations), 

ImVg,iz2, r, a) = -^^^i^ p^{b, z^) . (45) 

The evolution equation (|19D with the potential Vqq{z2,r2,a) containing this imaginary part 



was used in ||5^ , pO| , and nuclear shadowing in deep-inelastic scattering was calculated in good 
agreement with data. 

The analytical solution of Eq. ( pi| ) is only known for the harmonic oscillator potential V{r) oc 
r^. To keep the calculations reasonably simple we are forced to use the dipole approximation. 



(^qq{r,s) = C{s)r^' 



(46) 



The energy dependent factor C{s) is adjusted to reproduce correctly nuclear effects in the limit 
of very long CL Ic ^ Ra (the so called "frozen" approximation) , when 



Gqq{zi,ri; Z2, r2) => S{ri - r2) exp 



■2^99(^1) / dzpA{b,z) 



(47) 
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where the dependence of the Green function on impact parameter is dropped. The details are 
described in Appendix A. 

With the potential Eqs. p5| ) - p6| ) the solution of Eq. (|1^) has the same form as Eq. (|2T|), 



except that one should replace u ^ Q, where 



\/a'^(a) — i pAib, z) u a (1 — a) C(s) 

Q = -^-— ^ ^ -^ -—- . (48) 

u a{l — a) 

Guided by the uncertainty principle and the Lorentz transformation one can estimate the 
coherence time as in Eq. (^, where the effective mass of the qq pair is replaced by the vector 
meson mass. One can see the presence of a coherence length in the kinetic term of the evolution 
equation Eq. (|1^). Indeed, the effective mass squared of a qq pair is M|^ = (m^ + k^)/a{l — a). 
This is what the kinetic term consists of when the transverse momentum squared of the quark 
is replaced by k^ =^ A^. This dynamically varying effective mass controls the CL defined by 
the Green function, as compared to the oversimplified Eq. (|]) for the CL as given by the fixed 
mass mv- One can explicitly see the static part Q^ + ?Ti^/a(l — a) of the coherence length in 
the last phase shift factor in the Green function in Eq. (|21|). 

Depending on the value of Ic one can distinguish different regimes: 

(i) The CL is much shorter than the mean nucleon spacing in a nucleus {Ic — * 0). In this 
case G{z2,r2', Zi,ri) — ^ S{z2 — zi) since strong oscillations suppress propagation of the qq over 
longer distances. In this case the formation time of the meson wave function is very short as 
well, since it is described by the same Green function and is controlled by the formation time 
as given in Eq. (|^). Apparently, for light vector mesons // ~ Ic, so both must be short. In this 
case nuclear transparency is given by the simple formula Eq. (^) corresponding to the Glauber 
approximation^. 

(ii) In the intermediate case Ic -^ 0, but // ~ Ra, which can only be realized for heavy 

flavor quarkonia, the formation of the meson wave function is described by the Green function 

■^Note that the optical approximation is used throughout this paper only for the sake of easy reading. For 
numerical calculations we replace the exponential by a more realistic expression, exp(— cr Ta) => {1~(jTa/A). 
This has also been done in all of our previous publications, contrary to what is stated in ||55|. 
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and the numerator of the nuclear transparency ratio Eq. (^) has the form [Q , 



M 



7*A^VX 



>,Q^ 



lc^O;lfr^RA 



cfb dz pA{b, z) Fi{b, z) 



(49) 



where 



Fi{b,z) = I da I d^rid^r2'^'y{f2,a) G{z' ,f2; z,fi) (Tgg{ri, s) '$^*{fi,a) 



(50) 



This expression is illustrated in Fig. ^. The photon creates at the point z a colorless qq pair 





a b 

Figure 3: The incident virtual photon produces incoherently at the point 
z (quasielastic scattering) the colorless qq pair which then evolves prop- 
agating through the nucleus and forms the V-meson wave function (a). 
Alternatively, the photon can first produce diffractively and coherently 
at the point Z\ the colorless qq which then experiences quasielastic scat- 
tering at the point z (h). Propagation of the qq pair is described by the 
Green function (shaded areas). 



with transverse separation fi. The quark and antiquark then propagate through the nucleus 
along different trajectories and end up with a separation r2. The contributions from different 
paths are summed up giving rise to the Green function G{z',f2;z,fi) which is convoluted in 
(|4^) with the wave functions of 7* and V. This is the path integral technique suggested in 0. 
(iii) Ic ^ Ra (in fact, it is more correct to compare with the mean free path of the qq in a 
nuclear medium if the latter is shorter than the nuclear radius). In this case G{z2, r2\ zi, fi) — > 
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5{r2 — ri), i.e. all fluctuations of the transverse qq separation are "frozen" by Lorentz time 



dilation. Then, the numerator on the r.h.s. of Eq. (|4^) takes the form 



M 



7*A-^V'X 



s.Q' 



Ic:s>Ra 



d%TA{h) 



(fr / da 
Jo 



(51) 



X '$y{f,a)aqg{r,s) exp 



-2^mir,s)TA{b) 



^^.(f,a,Q^ 



In this case the qq attenuates with a constant absorption cross section like in the Glauber 
model, except that the whole exponential is averaged rather than just the cross section in 
the exponent. The difference between the results of the two prescriptions are the well known 
inelastic corrections of Gribov [0. 

(iv) The main and new results of the present paper address the general case with no 
restrictions for either l^ or If. No theoretical tool has been developed so far beyond the limits 
(i) - (iii) discussed above neither of which can be applied to electroproduction of light vector 
mesons at the medium high energies of HERMES and JLab. 

Even within the VDM the Glauber model expression interpolating between the limiting 
cases of low [(i), (ii)] and high [(iii)] energies has been derived only recently [0. We generalize 
that formalism to the LC dipole approach, and the incoherent photoproduction amplitude is 
represented as a sum of two terms |]SB| illustrated in Fig. ^ 



M 



■y*A-^VX 



:s,Q' 



d% / dzpAib,z) Fi{h,z)-F2{h,z 



(52) 



The first term Fi {b, z) introduced above in Eq. (|50D is represented by Fig. |^. Alone it would 
correspond to the short Ic limit (ii). The second term -F2(6, z) in ( ^2|) corresponds to the situation 
illustrated in Fig. ^. The incident photon produces a qq pair diffractively and coherently at the 
point Zi prior to incoherent quasielastic scattering at point z. The LC Green functions describe 
the evolution of the qq over the distance from Zi to z and further on, up to the formation of 
the meson wave function. Correspondingly, this term has the form. 



1 "" "^ 

F2{b,z) = ^ / dzipA{b,Zi) da d'^rid'^r2d'^r'^y{f2,a) 

-oo 

X G{z' -^ oo, f2; z, r) crgg{r, s) G{z, r; Zi, n) (Tgg{ri, s) ^^.(fi, a) 



(53) 
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Eq. (0) correctly reproduces the limits (i) - (iii). Indeed, at /c ^ the second term -^2(6, z) 
vanishes because of strong oscillations, and Eq. (|5^) reproduces the Glauber expression Eq. (§). 
On the other hand, at Ic ^ Ra the phase shift in the Green functions can be neglected and 
they acquire the simple form 0(22,^2; zi,fi) -^ 6{f2 — ri). In this case the integration over 
longitudinal coordinates in Eqs. ( |50D and (^) can be performed explicitly and the asymptotic 
expression Eq. ( pTD is recovered as well. Moreover, if one uses a constant dipole cross section 
^qq{.P) = ^Yot ^ then Eq. (0) recovers the general Glauber expression^ derived in |16|. 

3.2 Data for incoherent production: CT or coherence? 

Exclusive incoherent electroproduction of vector mesons off nuclei has been suggested in [^ 
as a sensitive way to detect CT. Increasing the photon virtuality Q^ one squeezes the produced 
qq wave packet. Such a small colorless system propagates through the nucleus with little 
attenuation, provided that the energy is sufficiently high {If ^ RjC) the fluctuations of the qq 
separation are frozen by Lorentz time dilation. Thus, a rise of nuclear transparency Tr^^'^{Q'^) 
with Q^ should signal CT^ Indeed, such a rise was observed in the E665 experiment at Fermilab 
for exclusive production of p° mesons off nuclei by a muon beam. This has been claimed in [31 
to be a manifestation of CT. 

However, one should be cautious to avoid mixing up the expected signal for CT with the 



effect of coherence length [63, 16 . Indeed, if the coherence length varies from long to short 



compared to the nuclear size the nuclear transparency rises because the length of the path in 



*Note that Eq. (|5^) and its Glauber model analog in |l6| include all coherent multiple scattering terms, 

contrary to a statement made in [ p5[ . 

^This process has a definite advantage compared to quasielastic electron scattering (e,e'p) suggested in [p^, p^ 

as a probe for CT. Indeed, in the latter case the energy of the photon correlates with its virtuality, v k, 2miyQ'^, 

and one has to increase Q^ just in order to increase i^ and keep the size of the ejectile "frozen" . This leads to a 

substantially diminished cross section, which is why no CT signal has been detected in this reaction so far (it 

is still possible to observe CT in this reaction at low energy studying the asymmetry of the quasielastic peak 

as function of xbj 16 0| ). In contrast, no correlation between v and Q^ exists in exclusive electroproduction of 

vector mesons. 



24 



nuclear matter becomes shorter and the vector meson (or qq) attenuates less. This happens 
when Q^ increases at fixed v. One should carefully disentangle these two phenomena. 

Long CL. It has been checked in |^3| that the coherence length at the kinematics of the 
E665 experiment is sufficiently long to neglect its variation with Q^ and to use the "frozen" ap- 
proximation, except at the highest values of Q^ > 5 GeV^. We calculated nuclear transparency, 
Tr"^^^, of incoherent (quasielastic) p^ production using Eq. (|52D and the simplified "frozen" ap- 
proximation Eqs. (1^) - (0). The results are depicted in Fig. || by solid and dashed curves 
respectively. One can see that fiuctuations of the size of the qq pair become important only 

1 



^A 



inc 



0.1 



— I 1 r- 



— I 1 r- 




n E665 



_] I I L. 



2 4 6 8 10 

.2 r r. ^ /2^ 



0.02 

2 4 6 

Q^ (GeV^ 

Figure 4: Q^- dependence of nuclear transparency for lead and calcium 
Trph and Tvca vs. The experimental points are from the E665 experi- 
ment [611. Both the curves and data for lead are rescaled by the factor 
1/2. Solid and dashed curves show our results using the LC Green func- 
tion approach Eq. ( ^ ) and the "frozen" approximation Eq. ( pl| ) respec- 
tively. 



at high Q^ causing a separation of the solid and dashed curves. At smaller Q^ the observed 
variation of Tr^^'^{Q'^) is a net manifestation of CT. The agreement with our model is surpris- 
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ingly good for calcium, while we underestimate the nuclear transparency at small Q^ for lead. 
This may be a manifestation of large Coulomb corrections as found in |^, which are of the 



order aemZ ~ 0.6 for lead. These corrections lead to a considerable deviation from the Born, 
one-photon approximation employed in [^] in order to obtain data for 7* A — > p° X (depicted 



in Fig. 1^) from raw data for /i A — i> ^' p^ X. This important problem needs further study. 

Medium long CL. The same process of incoherent electroproduction of p° is under study 
at lower energies, in the HERMES experiment at HERA and at JLab. In this case one should 
carefully discriminate between the effects of CT and CL [^, |16[. A simple prescription to 
eliminate the effect of CL from the data on the Q^ dependence of nuclear transparency is to 
bin the data in a way which keeps Ic = const. It means that one should vary simultaneously z/ 
and Q^ maintaining the CL Eq. @ constant, 

i^=llc{Q' + ml) . (54) 

In this case the Glauber model predicts a Q^ independent nuclear transparency, and any rise 
with Q2 would signal CT g. 

The LC Green function technique incorporates both the effects of coherence and formation. 
We performed calculations of Tr^^'^{Q'^) at fixed l^. starting from different minimal values of v, 
which correspond to real photoproduction in Eq. (|5^), 

The results for incoherent production of p and at Umin = 0.9, 2, 5 and 10 GeV {Ic = 
0.6—6.75 fni) are presented in Fig. ^jfor nitrogen, krypton and lead. We use the nonperturbative 
LC wave function of the photon with the parameters of the LC potential ao,i fixed in accordance 
with Eq. (RTf) at f = 1/2. The u and d quarks are assumed to be massless, but we use 



rris = 0.15 GeV. Nuclear transparency for (p is stronger than for p as one could expect, but the 
difference is not significant. In what follows we discuss only our results for p. 

For p-mesons the predicted variation of nuclear transparency with Q^ at fixed Ic is much 
stronger than was found in 0. Those calculations have been done in the hadronic representa- 
tion which is quite challenging due to the necessity to know all the diagonal and off diagonal 
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Figure 5: Q^ dependence of the nuclear transparency Tr^^'^ for exclusive 
electroproduction of p (left) and (right) mesons on nuclear targets ^^N, 
^^Kr and ^^^Pb (from top to bottom). The CL is fixed at U = 0.60, 1.35, 
3.37 and 6.75 fm. 
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diffractive amplitudes for the vector meson and its excitations, as well as all photoproduction 
amplitudes. The predictions made in f^ were based on the two-coupled-channel model without 
any estimate of the accuracy of such an approximation. According to quark-hadron duality 
the LC Green function method is equivalent to the exact solution of the general multi-channel 
problem in the hadronic representation. The comparison therefore demonstrates that the two- 
channel approximation substantially underestimates the effect of color transparency. 



To see the scale of the theoretical uncertainty of our model |2^ for nonperturbative effects 
we compare in Fig. ^ the results for the p-meson obtained using the nonperturbative (solid 
curves) and perturbative photon wave functions Eq. ([I^) with rrig = 0.15 GeV (dashed curves). 
The difference between the two sets of curves is insignificant. 
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Figure 6: The same as in Fig. |^ for lead, but calculated with both non- 
perturbative (solid curves) and perturbative (dashed) wave functions of 
the photon. 



Motivated by the too weak signal predicted for CT it was suggested in that instead one 
can study the effect of coherence which has never been observed experimentally. Indeed, it was 
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found in [34] that data for nuclear transparency for p production plotted as function of Ic well 
agree with what was predicted in |T^ to be the effect of the CL. Now we find a rather strong 
signal of CT which may also affect the Ic dependence of Tr'^^'^ and cause a deviation from the 
Glauber model expectations. We therefore revise the previous conclusions [^, |6^ . 

In the VDM-Glauber model nuclear transparency is a function of l^ only (neglecting the 
weak energy dependence of cr[of)^ however it becomes a function of two variables, Tr^^'^{lc, Q"^), 
as soon as CT effects are involved. Therefore, our current predictions for the Ic dependence 
of Tr^^'^ vary with Q^. They are plotted by dashed curves in Fig. |^ for different fixed values 
of Q^ = 0.5, 1, 2, 3, 5 GeV^ (from bottom to top) for nitrogen and krypton (left and right 
boxes, respectively). The nonperturbative wave function of the photon was used as for Fig. |. 
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Figure 7: Nuclear transparency for incoherent electroproduction of p 
off nuclei, nitrogen and krypton, as function of Ic at fixed Q^ = 
0.5, 1, 2, 3, 5 GeV^. Tlie solid curve is calculated at tlie mean val- 
ues of Ic and Q^ corresponding to each experimental point 



We do not show the results obtained with the perturbative wave function since they are pretty 
much the same, except in the region of small Q^ and short Ic where they are about 10% lower 
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that the nonperturbative results. 

The experimental points for nitrogen |Q which are plotted in Fig. correspond to different 



mean experimental values of Q^ [^]. This Q^ — lc correlation is incorporated in our calculations, 
and the results depicted by the solid curve agree well with the data. 

We thus arrive at the conclusion that the two quite different approaches, the VDM based 
Glauber model and QCD based LC Green function formalism, both provide good agreement 
with the HERMES data. This could not be possible if the data were plotted as function of Ic 
at fixed Q^. The observed agreement with the Glauber model seems to be accidental and a 
result of the Q^ — Ic correlation in the data. 

In order to discriminate between the two approaches one should plot the data differently. 
Fig. gives hope that the data are sufficiently accurate to detect a signal of CT if they are 
properly analyzed. Also additional data for krypton should soon become available from HER- 
MES. 

The expected signal for CT is a nonzero derivative d ln[TrA{Q'^)]/dQ'^ which is predicted 
in Fig. ^ to be similar for different nuclei and different values of Ic- One can make use of this 
fact and perform a common fit to all available data with only one parameter which is the slope 
of the Q^ dependence of nuclear transparency. The value of the logarithmic slope for the mid 
values, Q^ ~ 1 — 2 GeV^, of the HERMES kinematical range for p production is expected to 
vary within the interval 



dTr'riQ^'^ 0.07-0.11 GeV^^ for ^^N 



r-^ 



i,-const 0.14-0.17 GeV"^ for ^'^Kr 



(56) 



for Ic = 0.60 — 6.75 fm. Similar, but somewhat smaller values of the logarithmic Q^-slope are 
expected for 0. 

The curves in Fig. ^ demonstrate an interesting property. The slope of the Q^ dependence is 
steeper at small Q^ and Ic- For instance, the logarithmic derivative Eq. (^) equals 0.09 at Ic = 
0.6 fm, but is smaller, 0.07 at Ic = 1.35 fm This fact might be in variance with naive intuitive 
expectations. Indeed, Ic = 0.6 fm is short compared to th mean spacing of the bound nucleons. 
Since // ~ Ic at low Q^ one might expect the Glauber model to be a good approximation in this 
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case. Apparently, this is not the case, Fig. ^ demonstrates a steepest growth of Tr^^'^{Q'^) in 
this region. One can understand this as follows. If Ic is long, like in Fig. ^, then the formation 
length is long too. If >, Ic '^ Ra, and nuclear transparency rises with Q^ only because the 
mean transverse separation of the qq fluctuations decreases. If, however, Ic < Ra and fixed, the 
photon energy rises with Q^ according to Eq. (Q) and the formation length Eq. (|^) rises as 
well. Thus, these two effects, the Q^ dependence of // and the qq transverse size, add up and 
lead to a steeper growth of Tr^^^{Q^) for short Ic- 

One should conclude from this consideration that the CT effects are more pronounced at 
low than at high energies. This observation adds to the motivation for experimental searches 
for CT at HERMES and JLab. 

We also calculated the energy dependence of nuclear transparency at fixed Q^. The results 
for nitrogen and lead are shown by dashed curves in Fig. H for different values of Q^. The 
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Figure 8: Nuclear transparency for incoherent electroproduction '~f*A — > 
p^A as function of energy at Q^ = 0, 1, 3, 5, 10 GeV^ for nitrogen and 
lead. The solid and dashed curves correspond to calculations with and 
without gluons shadowing, respectively. 

interesting feature is the presence of a maximum of transparency at some energy. It results 
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from the interplay of coherence and formation effects. Indeed, the FL rises with energy leading 
to an increasing nuclear transparency. At some energy, however, the effect of CL which is 
shorter than the FL, is switched on leading to a growth of the path length of the qq in the 
nucleus, i.e. to a suppression of transparency. The maxima in the Ic dependence of nuclear 
transparency depicted in Fig. |^ are of the same nature. This also explains the unusual ordering 
of curves calculated for different values of Ic as is depicted in Fig. ^ 

4 Coherent production of vector mesons 

4.1 The formalism 

If electroproduction of a vector meson leaves the target intact the process is usually called 
coherent or elastic. The mesons produced at different longitudinal coordinates and impact 
parameters add up coherently. This fact considerably simplifies the expressions for the cross 
sections compared to the case of incoherent production. The integrated cross section has the 
form, 



^coh ^coh 



'^^*A-^VA 



2 

jd?q jd%e'^-'M''^Z\^yAh) = j dH\M';,\^yAh)W (57) 



where 

M'^,\^yAh) = J dzpAib, z) F,{b, z) , (58) 

— oo 

with the function Fi{h,z) defined in (^). 

One should not use Eq. (^^ for nuclear transparency any more since the t-slopes of the 
differential cross sections for nucleon and nuclear targets are different and do not cancel in 
the ratio. Therefore, the nuclear transparency also includes the slope parameter B^*n for the 
process '~^* N ^>V N , 

One can also define a t-dependent transparency for coherent electroproduction of vector 
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mesons, 



TrTit) 



da'X^/dt 



A^daN/dt\t=o 
where the differential cross section for coherent production 'y* A ^ V A reads 



(60) 



da"^'' 



dt 16 vr 



rf^fee*^^ / dz pA{b, z) Fi{b, z) 



(61) 



with Fi{b, z) defined in (|50D . This expression is simphfied in the hmit of long coherence time 
(t = V), 



rfaX^ 



dt 



1 
,:^Ra 4 vr 

1 



d%e'^-'^ I d^r{l - exp 

2 



- 2^5g(^s)^W 



X / cia\l'5^(f, a) \l'^.(r, a) 



(62) 



a form which resembles its VDM analogue ||TD|. 



4.2 Comparison with data and predictions for coherent production 



Using Eq. ( ^9]) we can also calculate the normalized ratio of coherent cross sections on two 
nuclei i?coh(^iM2) = Tr'Xl'/Tr''^^. The results of calculations for R^oh{Ph/C) and Rcoh{Ca/C) 
are depicted by solid curves in Fig. |^ as well as corresponding data from the E665 experiment 
|]6T| shown by squares and triangles, respectively. We performed calculations of Tr"^^ at mean 
photon energy z/ = 138 GeV with the Q^ dependent slope given by Eq. (^Tj). All effects of CL 
and CT are included via the LC Green function formalism. For such a high energy one can think 
that the "frozen" approximation Ic ^ Ra is good. In order to check how variation of the CL 
affects the nuclear transparency we repeated our calculations in the "frozen" approximation and 
plotted the results as dashed curves in Fig. |^. We see that the accuracy of this approximation 
is rather good for calcium, while for lead it significantly deviates from the exact result at 
Q^ > 2 GeV^. The reason is obvious, the heavier the nucleus, the less the approximation 
Ic ^ Ra is fulfilled. We also see that the contraction of the CL with Q^ causes an effect 
opposite to CT, namely nuclear transparency is suppressed rather than enhanced. Therefore, 
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Figure 9: Q^- dependence of the total cross section ratio Rcohi^/C) = 
12(j'^^ /Aa^^ for the coherent process 7*A —^ p^A. Experimental points 



are from E665 j|6l|/ for Ph/C (squares) and Ca/C (triangles). Solid curves 
include the variation of Ic and If with Q^. Dashed curves are calculated 
in the "frozen" approximation l^ ^ Ra- 



there is no danger that CL effects can mock CT, and one may think that this is an advantage of 
coherent compared to incoherent production ||6^. However, at medium energy the suppression 
of nuclear transparency at short CL is so strong that no rise of nuclear transparency with Q^ 
might be observable. 

Note that in contrast to incoherent production where nuclear transparency is expected to 
saturate as Tr™'^((5^) ^ 1 at large Q^, for the coherent process nuclear transparency reaches 
a higher limit, Tr'^^(Q'^) -^ A^^^ (of course, A^^^ is valid only for very large nuclei, otherwise 
it is an approximate number). The dashed curves in Fig. ^ nearly reach this upper limit at 
g2 ^ 10 GeV^. 

One can eliminate the effects of CL and single out the net CT effect in a way similar to 
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what was suggested for incoherent reactions by selecting experimental events with l^ = const. 
We calculated nuclear transparency for the coherent reaction '-f*A —>■ p{(f))A at fixed values of 



Ic- The results for Ic = 1.35, 3.37, 13.50 fm are depicted in Fig. M for several nuclei. The effect 
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Figure 10: Tie same as in Fig. |^, but for coherent production of p and (j), 
YA -> VA. 
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is sufficiently large to be observable, the logarithmic derivative varies within the interval, 



1 



Tr'^^iQ'^ 



dTr''X^{Q'^) 0.14-0.07 GeV"2 for ^^A^ 

lc=const ^ ^ ^ n 1 fr /^ „■» r-2 f„„ 84 i 



dQ^ 



(63) 



0.10-0.15 GeV"^ for ^^Kr 



for l(. = 1.35 — 13.5 fm. Again, like in the case of incoherent production, the logarithmic 
derivative decreases at large Ic- The magnitude of the expected CT effect is similar to the value 
predicted for the incoherent production in Eq. (|56[) and is slightly smaller for phi than for p. 

We also calculated nuclear transparency as function of energy at fixed Q^. The results 
for p produced coherently off nitrogen and lead are depicted by dashed curves in Fig. |TT] at 
Q^ = 0, 3, 10 GeV^. Tr'^'^ is very small at low energy, what of course does not mean that 
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Figure 11: Nuclear transparency for incoherent electroproduction 7*A — ;> 
pA as function of energy at Q^ = 0, 3, 10 GeV^ for nitrogen and lead. 
Tie solid and daslied curves correspond to calculations with and without 
gluon shadowing, respectively. 

nuclear matter is not transparent, but the nuclear coherent cross section is suppressed by the 
nuclear form factor. Indeed, the longitudinal momentum transfer which is equal to the inverse 
CL, is large when the CL is short. However, at high energy Ic ^ Ra and nuclear transparency 
nearly saturates (it decreases with u only due to the rising dipole cross section). The saturation 
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level is higher at larger Q'^ which is a manifestation of CT. 



4.3 Transverse momentum distribution 



Another manifestation of CT is a modification of the diffractive pattern in the momentum 



transfer dependence of the coherent cross section ||55|. Indeed, the effect of CT on the nu- 
clear transparency depends on impact parameter. Fourier transformation of such a modified 
amplitude will apparently result in a shifted positions of the diffractive minima. Indeed, calcu- 
lations performed in the "frozen" approximation assuming sufficiently high energy (/c ^ Ra) 
lead to the t-dependence of nuclear transparency from Eq. (|60D depicted in Fig. |12|. We see 
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Figure 12: Nuclear transparency for coherent electroproduction 7*^4 — >■ 
p^A as function of momentum transfer squared calculated for nitrogen 
and lead in the limit oflc ^ Ra- The solid and dashed curves correspond 
to Q^ = and 10 GeV^ respectively. 



that the CT effects shift the position of the diffractive minima to larger t. To understand the 
sign of the effect, we can use the approximate dipole cross section crgg(r) = Cr'^. Further, we 
can approximate the product of the photon and vector meson wave functions by a Gaussian 
oc exp(— r^/(r^)). The partial amplitude (amplitude for given impact parameter) of elastic 



37 



production 7* -^ p° takes the form , 

C{r')TAib) 



f (fr ^*y{r} ^^, (f , g2) j 1 _ exp 



-Icr'T^ib) 



(64) 



2 + C(r2)T^(6) 

The mean size (r^) of the qq wave packet decreases with Q^ suppressing the partial amphtude 
from Eq. (|64D. It follows from (0) the suppression is smaller on the periphery of the nucleus 
than in the center. This implies that the slope of the t-distribution B^*a = (^^)/2 should 
decrease with Q"^, i.e. the minima should move to larger values of t. 

Although such a modification of the diffractive pattern should signal CT, the effect is very 
weak and its observation does not look feasible. 

5 Other pitfalls in the search for CT at low energies 

In order to avoid the effect of CL which leads for incoherent production to a nuclear transparency 
rising with Q"^ and mimics CT we suggested in Sect. ^^ to study the Q^ dependence in data 
samples which are preselected to have the same Ic- As soon as the CL effects are excluded, 
the Glauber model predicts no variation of the transparency with Q^. There are, however, still 
other effects, not related to CT, which cause a growth of CT even if /^ = const. 

5.1 Standard inelastic corrections 

The effect of CT can be treated in the hadronic representation as a multichannel problem 
10, |6^, 1^: the incident virtual photon produces diffractively on a bound nucleon either the 



ground state V, or any excitation. Only this stage of the process is Q^ dependent. The produced 
states propagate further through the nucleus experiencing diagonal and off diagonal diffractive 
transitions. Eventually, the state V is detected at macroscopic distances. These modifications 
of the Glauber single-channel approximation are at the heart of Gribov's inelastic shadowing. 
The miracle of CT is the expectation that all those large amplitudes must cancel leaving 
only one amplitude, namely the direct production of the V. There is no hint from the hadronic 
representation that this should happen. We have no data for most of those amplitudes and no 
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hope to measure them in future. Only the gauge invariance of QCD dictates this very nontrivial 
behavior which is not present in any of the old fashion models (colorless constituent quarks, 
etc.). 

Nevertheless, we do have data for single diffraction which allow to calculate some of the 
lowest order inelastic corrections. Although these corrections are part of the whole CT phe- 
nomenon they are model independent (provided that those models are fitted to available data). 
In particular, the nuclear medium is known to be more transparent than expected using the 
Glauber model [^. Indeed, if the produced V state experiences inelastic diffraction inside the 
nucleus, it is gone from the detected channel according to the rules of the Glauber approxi- 
mation. However, there is still a possibility to recover and come back to the V channel in a 
subsequent collision, as is illustrated in Fig. |l^. Apparently, this process increases the survival 




Figure 13: Intermediate diffractive excitation of a vector meson produced 
by the incident virtual photon and propagating through the nucleus. 

probability for the V state. There is clear experimental evidence that this takes place. The to- 
tal cross sections of hadron- (neutrons, neutral kaons) nucleus interactions measured with high 
accuracy [pTl pSf are smaller than the Glauber model predictions. This deviation increases 



with energy as it is controlled by the nuclear form factor which depends on the longitudinal 
momentum transfer qi in the diffraction dissociation. 

Therefore, transparency of nuclear matter for hadrons increases with energy and this fact 
leads to a rising Q^ dependence if data are selected according to the condition Ic = const. 
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Indeed, the energy rises according to the correlation v = [Q^ + my)/2lc. Of course this effect 
is well known for the total cross sections since the pre-QCD era. It cannot be (and never was) 
interpreted as a manifestation of CT. Although these inelastic corrections are part of the CT 
phenomenon, one should admit that they exist independently of the answer to the question 
whether the CT is true or not. 

Thus, one should be cautious in interpreting a rise of nuclear transparency as function of Q^ 



at fixed Ic- The correction under discussion calculated in ||69| for quasielastic high-p^ electron 
scattering, A{e, e'p)X was found to be indistinguishable from the predicted CT effect up to 
rather high Q^ of a few tens of GeV^. 

The deviation of the transparency from the Glauber model prediction is calculated as ||69(] , 



inc//^2\ 



Tr^riQ 



(fb / dz pAib, z) exp 



a. 



VN 



dz' pA{b,z') 



X 



1 + An dM^ - ^ ' 



dW^dt 



FA{b,z,qL] 



t=o 



(65) 



Here F^ib, z, qi) in ( |5BD is the so called longitudinal form factor of the nucleus calculated at a 
given impact parameter b and production coordinate z, 



FA{b,z,qL) = / dz' pA{b,z') cos{qLz') . 



(66) 



We use the same parameterization for the single diffraction cross section as in [|67| , |69| except 
for the normalization which is reduced by the factor 2/3 as is suggested by the triple Regge 
phenomenology. Although it is a rather rough estimate, it is sufficient for our purpose, since 
the effect turns out to be very weak. The results of our calculations for the Q^ dependence 



of the transparency are depicted by solid curves in Fig. |Tj at fixed Ic = 1.35 and 6.75 fm (the 
bottom and upper curves respectively). One can see that although Tr^^^ grows with Q^, this 
effect when compared with Fig. ^ is too weak to be mixed up with the signal of CT. Indeed, the 
derivative d \n{Tr^]^'^)/dQ'^ evaluated at Q^ = 1 — 2 GeV^ equals 0.011 for nitrogen, 0.025 for 
krypton and 0.033 for lead. This is nearly an order of magnitude less than what was estimated 
in (|56D as a signal for CT. 
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Figure 14: The same as in Fig. ^ but without any CT effects. Only 
the standard inelastic corrections which make the nuclear matter more 
transparent are included. The dashed curves correspond to the Glauber 
approximation. Each couple of solid curves correspond to Ic = 1-35 fm 
(bottom curve) and Ic = 6.75 fm (upper curve). The two sets of curves 
correspond to nitrogen (top) and lead (bottom) 

5.2 Finiteness of the p meson lifetime 

Some of the effects have been calculated above at rather low energies when the lifetime of 
the p is comparable with the nuclear size. For instance, at z/ = 2 GeV the mean path length 
up to the decay is only 2.7 fm. Two pions have a smaller survival probability than the p 
therefore the nuclear transparency should be smaller than is expected within the Glauber 
approximation disregarding decays. However, as function of energy the decay path length 
increases and eventually the nuclear transparency must reach the value corresponding to the 
Glauber model. Again, at l^ = const the energy of p rises with Q^ and the nuclear transparency 
must grow with the decay length. This effect might cause a problem in identifying the signal of 
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CT at low energy, and it should be considered with care. This correction is of less importance 
for (j) production. 

We corrected the Glauber formula for the finite decay length in the following way, 

oo oo 

TrT = 7{iy) Jd% J dzipA{b,z)Jdz2 

— oo Zi 

Z2 

pN 



dz'pA{b,z') - 2 a^^ dz pA{b, z') + 'y{u){z2 - zi) 



(67) 



7(-) = -^^^ , (6^ 



^2 _ j^2 



X exp 
where 



p 
is the Lorentz enhanced decay length of the p meson, Tp = 0.15 GeV is the total decay width 

of the p. Like previously, u correlates with Q^ via Eq. (l5^). 

The results of our calculations are shown in Fig. |l^ as function of Q^ at different fixed Ic- 
The effect turns out to be very weak compared to the expected effect of CT demonstrated in 
Fig. ^. For a CL of Ic = 1.35 fm (6.75 fm) the results are shown by the bottom (upper) solid 
curves for both nuclei. The derivative d ln{Tr^^^)/dQ'^ evaluated at Q^ = 1 — 2 GeV^ equals 
0.014 (0.003) and 0.03 (0.008) for nitrogen and lead respectively. We conclude that the effect 
of the finite p decay length cannot be mixed up with a CT signal. 

It is disputable whether the two pions emerging from a p decay immediately starts to 
attenuate with twice the pion absorption cross section. One might think about two pions 
which strongly overlap at the production point, then their cross section should be reduced due 
to color screening. This is a general problem of how one should decide whether decay has 
already happened or not. It is easier to understand for the example of photon radiation by an 
electron. One can treat the photon as being originally a part of the static Coulomb field of the 
electron which is then shaken off by the interaction with a target. Only when the photon and 
electron become incoherent they start acting as independent partons. It takes, however, a time 
span dictated by the uncertainty principle to discriminate between a coherent system, electron 
and its field, and an incoherent pair of an electron plus a photon. This time is called radiation 
or coherence time. In analogy, one can say that the p meson has already decayed when the 
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Figure 15: The same as in Figs. ^ and |I^ but with no effects of either 
CT or any inelastic correction included. Only the finite p meson decay 
length which leads to an enhanced absorption is taken into account. For 
each nucleus the dashed curve corresponds to the Glauber approximation 
(Tp -^ 0), the bottom and upper solid curves correspond to Ic = 1.35 and 
6.75 fm respectively. 



two pions become incoherent. In this case they interact independently and Eq. (pTf) is vahd. 
However, while the pions are still coherent, one should treat them as intrinsic components of 
the p meson. 



6 Gluon shadowing 

At very small xbj the density of gluons should eventually deviate from the predicted by the 
DGLAP evolution (as we mentioned, whether an indication of this saturation effect was already 
seen at HERA is controversial). In the infinite momentum frame this phenomenon looks like 
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gluon-gluon fusion. It corresponds to a nonlinear term in the evolution equation |72|. This 
effect should lead to a suppression of small-xsj gluons in a nucleon. 

One may expect a precocious onset of the saturation effects for heavy nuclei. In the infinite 
momentum frame of the nucleus the gluon clouds of nucleons which have the same impact 
parameter overlap at small xbj in longitudinal direction. This allows gluons originated from 
different nucleons to fuse leading to a gluon density which is not proportional to the density of 
nucleons any more. This is gluon shadowing. 

Such a parton model interpretation is not Lorentz invariant, the same phenomenon looks 
quite different in the rest frame of the nucleus. It corresponds to the process of gluon radiation 
and shadowing corrections related to multiple interactions of the radiated gluons in the nuclear 



medium [|1^]. This is a coherence phenomenon known as the Landau- Pomeranchuk effect, 
namely the suppression of bremsstrahlung by interference of radiation from different scattering 
centers. It demands a sufficiently long coherence time of radiation, a condition equivalent to 
demanding a small Bjorken xbj in the parton model. 

Although the two interpretations look so different, one can get a hint that they are the same 
phenomenon relating them to the Reggeon graphs which are Lorentz invariant. The double- 
scattering correction to the cross section of gluon radiation depicted in Fig. p!6| a corresponds to 
the absorptive part of elastic pA amplitude shown in Fig. p!6|b. Since the initial and final nucleons 
are colorless, each pair of exchanged gluons attached to the same nucleon is in a colorless state. 



i.e. represents the Pomeron (in the Born approximation). Thus, the Feynman graph in Fig. p!q b 
is a part of the triple-Pomeron diagram shown in Fig. |TB|c. It can be interpreted as fusion of 
two Pomerons originated from different nucleons, 2IP -^ IP. This observation bridges the two 
interpretations of gluon shadowing. 

Note that in the hadronic representation such a suppression of parton density corresponds 
to Gribov's inelastic shadowing |l[] which is related to the single diffraction cross section. In par- 
ticular, gluon shadowing corresponds to the triple-Pomeron term in the diffractive dissociation 
cross section which enters the calculations of inelastic corrections (Fig. |l6|c). 

There are still very few numerical evaluations of gluon shadowing in the literature, all done 
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a b c 

Figure 16: The double scattering correction to gluon radiation in the rest 
frame of the target nucleus (a) . The absorptive part of the corresponding 
elastic pA amplitude (b). The triple-Pomeron graph representing fusion 
2IP ^IP (c). 



in the rest frame of the nucleus using the idea from [|T9l . It turns out that even at low Q^ in 
the nonperturbative domain where one should expect the strongest shadowing effects they are 
rather weak. Indeed, data for diffractive excitation of the incident hadrons to the states of large 
mass, the so called triple-Pomeron region, show that the cross section is amazingly small, an 
order of magnitude smaller than one could expect by comparing with the cross section of small 



mass excitation |T8|. To explain such a smallness one has to assume a rather small radius of 



propagation of the LC gluons, tq ~ 0.3 fm |g3|, |7l[]. It is clear that such a small quark-gluon 
fluctuation also leads to a rather weak gluon shadowing. 

To incorporate the smallness of the size of quark-gluon fluctuations into the LC dipole 
approach a nonperturbative LC potential describing the quark-gluon interaction was introduced 
into the Schrodinger equation for the LC Green function describing the propagation of a quark- 
gluon system. The strength of the potential was fixed by data on high mass (Mx) diffraction 
pp —>■ pX |]23[. This approach allows to extend the methods of pQCD to the region of small Q^. 



Since a new semihard scale I/tq ~ 0.65 GeV is introduced one should not expect a substantial 
variation of gluon shadowing at Q^ < 4/rg. Indeed, calculations performed in p3| for Q"^ = 
and 4 GeV^ using different techniques led to about the same gluon shadowing. At higher Q^ 
shadowing slowly (logarithmically) decreases in accordance with expectations based on the 
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evolution equation W^ . 



We repeated the calculations [^ of the ratio of the gluon densities in nuclei and nucleon, 

^^^"-- ^ ) - AG,{xs,,Q^) ^ ' - -^zr ' ^ ^ 

where Aa{qqG) is the inelastic correction to the total cross section cr/ot related to the creation 
of a qqG intermediate state, 

Acr(ggG) = Re / dz2 / dzi pAib, zi) p^ib, Z2) I <f'X2<f'y2d^XicP'yi da^ 



Q 

-00 —00 



X F}^._g^(.{x2, m, aq, ac) GqqG{x2, m, Z2; xi,yi, zi) F^^^gqcixi, yi, aq, ac) ■ (70) 

Here x and y are the transverse distances from the gluon to the quark and antiquark, respec- 
tively, aq is the fraction of the LC momentum of the qq carried by the quark, and ac is the 
fraction of the photon momentum carried by the gluon. F^*^qqG is the amplitude of diffractive 
qqG production in a 7*A^ interaction p3|] . 



9 

F^,^qqG{x,y,aq,aG) = -'^qq{aq,x-y) 



(Tqq{x) + aqq{y) - aqq{x - y) 



(71) 



where ^qq and ^qc are the LC distribution functions of the qq fluctuations of a photon and qG 
fluctuations of a quark, respectively. 

GqqG{x2, y2, Z2] Xi, ^i, Zi) is the LC Green function which describes propagation of the qqG 
system from the initial state with longitudinal and transverse coordinates zi and Xi,yi, respec- 
tively, to the final coordinates (-22, 3^2,^2) • For the calculation of gluon shadowing one should 
suppress the intrinsic qq separation, i.e. to assume x = y. In this case the Green function 
essentially simplifies and describes propagation of a gluon-gluon dipole through a medium. 

An important finding of Ref. [^ is the presence of a strong nonperturbative interaction 
which squeezes the gluon-gluon wave packet and substantially diminishes gluon shadowing. The 
smallness of the gluon-gluon transverse separation is not a model assumption, but is dictated 
by data for hadronic diffraction to large masses (triple- Pomeron regime) which is controlled by 
diffractive gluon radiation. 
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Further calculational details can be found in [23]. Here we calculate Rq [Eq. (^)] for 
different nuclear thicknesses TA{b). Since we use an approximation of constant nuclear density 



(see [Appendix 7\| ) , TA^b) = po L, where L = 2 J R\ — &^, the ratio Rg{xbj, Q^) is also implicitly 
a function of L. An example for the calculated L-dependence of Rg{xbj, Q"^) at Q^ = 4 GeV^ is 
depicted in Fig. |1^ for different values of xbj- As one should expect, the longer L, the stronger 
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Figure 17: The ratio of nucleus to nucleon gluon densities as function of 
the thickness of the nucleus, L = T{h)/ pq, at Q"^ = A and different fixed 
values of Xbj- 
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is gluon shadowing at small xbj- 

We calculated the gluon shadowing only for the lowest Fock component containing just one 
LC gluon. In terms of the parton model it reproduces the effects of fusion of many gluons to 
one gluon (in terms of Regge approach it corresponds to nIP -^ IP vertex). Inclusion of higher 
multigluon Fock components is still a challenge. However, their effect can be essentially taken 
into account by eikonalization of the calculated Rg{xbj, Q"^)-, as argued in |73|. In other words, 
the dipole cross section which is proportional to the gluon density at small separations, should 
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be renormalized everywhere, 

(^qq => Rg CTqq ■ (72) 

Such a procedure makes the nuclear medium more transparent. This could be expected since 
Gribov's inelastic shadowing is known to suppress the total hadron-nucleus cross sections, i.e. 
to make nuclei more transparent [0, |69[] . 

It is interesting, that the cross section of incoherent electroproduction of vector mesons is 
rather insensitive to gluon shadowing. Indeed, although the renormalization Eq. (^) suppresses 
the pre-exponential factor aqq{r,s) on the r.h.s. of Eq. (pTD, it simultaneously increases the 
exponential. These two effects essentially cancel. Indeed, our predictions for the effect of gluon 
shadowing for incoherent production '~f*A -^ p^X depicted in Fig. ^ demonstrate a rather small 
difference between the curves with (solid) and without (dashed) gluon shadowing. 

A few observations are in order. First, the onset of gluon shadowing happens at rather 
high energy u > 100 GeV. This corresponds to the claim made in ||2^ that the onset of 



gluon shadowing requires smaller xbj than the onset of quark shadowing. This is because the 
fluctuations containing gluons are in general heavier than the qq and have a shorter CL. 

Then, one can see that a stronger effect of gluon shadowing is expected for nitrogen than 
for lead. Although it contradicts simple intuition, it is easy interpreted. The renormalization 
of the dipole cross section Eq. ([72| ) may either suppress or enlarge the incoherent cross section 
in Eq. (|5T| ) depending on the value of the nuclear thickness function T{b). Namely, it should 
lead to a suppression for small T{b), but to an enhancement for large T{b). Indeed we observe 
this trend in Fig. p. Some enhancement (antishadowing) can be seen for lead at Q^ = 0. The 
results presented for nitrogen show that the effect is maximal at intermediate values of Q^ 
while it is smaller at Q^ = and 10 GeV. This is a result of the same interplay between the 
pre-exponential factor and the exponential in Eq. (0). 

The implication of gluon shadowing for the case of coherent production 'y*A — > VA is 
clearer. It is easy to understand that it always diminishes the coherent cross section. In terms 
of VDM the photoproduction cross section is related to the total VA cross section which is 
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always reduced by inelastic corrections. One can also see from Eq. (^) that the suppression 
of the dipole cross section by the renormalization [Eq. (|72D ] can result only in a reduction of 
the cross section. Thus, we expect much stronger effects of gluon shadowing for coherent than 
for incoherent production. The results of calculations based on the exact expression Eq. (pl|) 
confirm this conjecture. The predicted suppression of the coherent cross section is stronger for 
heavy nuclei and low Q^ . 

7 Summary and conclusions 

Electroproduction of vector mesons off nuclei is subject to an interplay between coherence 
(shadowing) and formation (color transparency) effects. Conventionally, one can associate 
those effects with the initial and final state interaction respectively. We developed a rigorous 
quantum-mechanical approach based on the light-cone QCD Green function formalism which 
naturally incorporates these interference effects. Our main results and observations are the 
following. 

• The suggested approach allows to find for the first time a comprehensive model for the long 
standing problem of exclusive electroproduction of vector mesons off nuclei. The main 
result of the paper, Eq. (|52D, interpolates between the previously known low and high 
energy limits for incoherent production. Eq. ( ^8]) does the same for coherent production. 

• The onset of coherence effects (shadowing) can mimic the expected signal of CT in inco- 
herent electroproduction of vector mesons at medium and low energies. In order to single 
out the formation effect data must be binned in Ic and Q^. Observation of a rising nuclear 
transparency as function of Q^ for fixed Ic would signal color transparency. 

• Due to quark-hadron duality the Green function formalism under consideration is equiv- 
alent to a solution of the full multi-channel problem in the hadronic representation. We 
found a much steeper Q"^ dependence of nuclear transparency (Fig. |[), i.e. a stronger 
signal of CT, than was predicted in [^ within the two-channel approximation. Moreover, 
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the slope of the Q^ dependence is even higher at lower energies. This should allow to 
detect a signal of CT in experiments with the HERMES spectrometer and especially at 
JLab. 

The successful experimental confirmation [Q of the predicted coherence length effects 



1 16] seems now to be a accidental consequence of the specific correlation between Q^ and Ic 
in the HERMES data. The present parameter-free calculations well describe the observed 
variation of nuclear transparency with Ic (Fig. |^) as a result of a complicated interplay 
between the effects of CT and coherence length. 

There are other effects which may cause a rise of nuclear transparency with Q^ at l^. = 
const thus mimicking a signal of CT. These are the lowest order inelastic corrections 
which are rather precisely fixed by available data for diffraction, and the finite lifetime 
of vector mesons. Both effects lead to the more transparent nuclear medium at higher 



energies, i.e. at larger Q due to the correlation in Eq. (|5^) between u and Q . We found 



however, both effects to be too weak (Figs. [T^ and |T5D to be relevant 



The effects of CT in coherent production of vector mesons are found to be less pronounced. 
Although transparency decreases with Q^ and does not mimic CT in this case, the CL 
effects significantly modify the Q^ dependence and may completely eliminate any signal 
of CT at medium energies. Besides, the cross section of coherent production is very small 
at low energies (Fig. |TI|). 

The effects of CT modify the impact-parameter dependence of the amplitude of coherent 
production by diminishing the mean square of the impact parameter in the interaction 
amplitude. Therefore the positions of the diffractive minima in the differential cross 



section are expected to shift to larger values of |t| (Fig. ^^. However the effects seems to 
be too small to be reliably observed. 

Although it has been known how to calculate nuclear transparency in the high energy 



limit Ic ^ Ra m, 0]' ^^^ effect of gluon shadowing was missed. We calculated nuclear 
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suppression of gluons at small xbj within the same LC approach treating it as shadowing 
corrections for the higher Fock states containing gluons. The nonperturbative interac- 
tion of the LC gluons significantly reduces the predicted magnitude of gluon shadowing 
(Fig. |1^. Although the amplitude of meson production off a bound nucleon is suppressed 



due to a reduced amount of gluons in the nucleus, the same effect makes the nuclear 
medium more transparent and enhances the meson survival probability. For incoherent 
p production these two effects nearly compensate each other for heavy nuclei (Fig. H). 
The cross section for coherent production is less for more transparent nuclei, therefore 
the effect of gluon shadowing is more pronounced (Fig. |Tl|). These corrections are not 
important at HERMES or JLab energies, but are significant at the higher energies of 
eRHIC and for coherent Coulomb production in heavy ion collisions at RHIC. 

Concluding, the predicted large effects of CT in incoherent electroproduction of vector 
mesons off nuclei open new possibilities for the search for CT with medium energy electrons. 
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Appendix A Approximation for the dipole cross section 

To keep the calculations simple we are forced to use the approximate dipole cross section 



Eq. (|U) which allows to obtain the Green function in an analytical form as is described in 
Sect. p.l| . We fix the factor C{s) by demanding that calculations employing the approximation 
Eq. (^61) reproduce correctly the results based on the realistic cross section in the limit Ic ^ Ra 
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when the Green function takes the simple form (47). Thus, for incoherent production of vector 
mesons the factor C{s) is fixed by the relation, 

2 



Jd%TA{b) 



J d r r exp 



-^CT,L{s)r'TA{b) 



T,L, 



S da-^*y{f,a)-^ I {f,a 



J d^r r2 / da ^^(f, a) ^'r{r, a 



Jd%TA{b) J d'^r agg{r, s) exp -\(Tqq{r,s)TA{h) 



^ da^!*y{r,a)^! r{f,a 



TXi 



J dh- aqq{r, s) J da ^|.(f, a) ^'r (f, a 



Correspondingly, for coherent production the factor C{s) is fixed by the relation. 



(A.l) 



Jd'b 



T,L, 



J d^r J da ^;^(f, a) %*{f, a) <1 - exp 



-lCTAs)r^TA{b) 



Jd% 



J d^r J da ^*y{f, a) Ct,l{s) ^^ ^r (^' « 
/ d^r J da ^^(r, a) *^i'^(r, a) 1 1 - exp -^ aqq{r, s) TA{b) 



.T,L, 



(A.2) 



/ dV / da ^*y{r, a) ^^^(r, s) %r{f, a) 

To take advantage of the analytical form of the Green function which is known only for the 
LC potential Eq. (^) with a constant nuclear density, we use the approximation pA{b,z) = 
Po 0(-Ra — b"^ — z"^)- Therefore we have to use this form for Eqs. ( [A.l|) and ( [A.2|) as well. The 
value of the mean nuclear density po has been determined using the relation. 



d^b 



1 - expl -cxq po J R\ - b^ 



d% 



1 — exp[ T(h) 



(A.3) 



where the nuclear thickness function TA{b) is calculated with the realistic Wood-Saxon form of 
the nuclear density. The value of po turns out to be practically independent of the cross section 
(To in the range from 1 to 50 mb. 
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